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Abstract 

We study the large N saddle points of weakly coupled M = 4 super Yang- 
Mills theory on S 1 x S 3 that are described by a commuting matrix model 
for the seven scalar fields {Aq,^j}. We show that at temperatures below 
the Hagedorn/'deconfinement' transition the joint eigenvalue distribution is 
S 1 x S 5 . At high temperatures T S> 1/R S 3, the eigenvalues form an ellipsoid 
with topology S 6 . We show how the deconfinement transition realises the 
topology change S 1 x S 5 — > S e . Furthermore, we find compelling evidence 
that when the temperature is increased to T = l/(\/ r \R s i}) the saddle with 
topology changes continuously to one with S 5 topology in a new second order 
quantum phase transition occurring in these saddles. 
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1 Introduction and summary 



In this work we will combine two strands of research that have each generated inter- 
esting results within the AdS/CFT correspondence [1] over the past couple of years. 

On the one hand, it has been understood how the dynamics of certain supersym- 
metric sectors of M = 4 super Yang-Mills theory on a spatial S 3 are described by 
matrix models for the scalar fields $j of the theory. Furthermore, it was discovered 
that the eigenvalue distributions of these scalars directly reconstruct a dual spacetime 
geometry [2, 3, 4, 5]. These works have provided a striking realisation of emergent 
geometry within the AdS/CFT setup. 

On the other hand, in the same theory at finite temperature it was shown that at 
weak 't Hooft coupling there is a phase transition at a critical temperature [6, 7]. This 
transition appears to be similar to the Hawking-Page transition at strong coupling, 
which describes the appearance of a black hole in the dual geometry [8, 9]. In partic- 
ular, for both phase transitions the eigenvalue distribution of the time component of 
the gauge field, A , becomes a nonuniform distribution on S 1 at the transition and 
ultimately becomes a gapped distribution as the temperature is increased. 

Combining the insights of these works, one might hope that studying the joint 
eigenvalue distribution of {A , $j} could lead to an understanding of the dual black 
hole geometry in the weak coupling regime. In fact, a fundamental question is whether 
or not there is a well defined sense in which it is useful to conceive of a weakly coupled 
plasma as being dual to a black hole. Various recent works have posed this as a 
question about correlators in real time physics [10, 11, 12, 13, 14]. The question is 
not purely of conceptual interest, as some discussion of the fireball created in the 
Relativistic Heavy Ion Collider attests, see for instance [15, 16]. The approach here 
will be Euclidean and hence concerned with equilibrium physics. 

In most works on the finite temperature theory at weak coupling, the scalar fields 
are integrated out on the grounds that they acquire thermal masses, as well as hav- 
ing a classical mass due to being conformally coupled to a spatial S 3 . However, at 
temperatures TR ~ A _1//2 , the one loop thermal mass of A is of the same magnitude 
as that for the scalar fields, so A and <3>j (and indeed the Ai [17]) should participate 
equally in the dynamics. The full dynamics of the theory is difficult to study, as 
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there are infinitely many coupled modes that can condense, and furthermore generic 
condensates will not commute. In this work we will make the drastic (but consis- 
tent) simplification of only considering saddles in which the homogeneous modes of 
{A Q , $j} condense and where the condensates commute. We believe that these are 
interesting saddles to consider; they are the most 'geometric' of possible saddles at 
weak coupling and can be described very explicitly. They may also have a connec- 
tion with the dominant geometrical saddles at strong coupling. However, it is clear 
[17] that they are not the absolute minima of the theory. Phase transitions in these 
saddles are not transitions of the full theory. 

The effective potential necessary for this study, the potential for commuting scalars 
and Aq, was computed in [18], who used it to show that the potential reveals a weak 
coupling analogue of the Gregory-Laflamme instability of the small black hole towards 
localisation on S 5 [18]. In [19] it was found that allowing one scalar field to be non- 
zero leads to a non-trivial joint eigenvalue distribution of {A ,Qj}. This is due to 
a logarithmic repulsion between the scalar eigenvalues at large N overcoming the 
classical and thermal mass terms. However, the solutions found in [19] break R 
symmetry in picking out a given scalar field. 

In this work we find the eigenvalue distribution of {Aq, $ j} that minimises the 
one loop effective potential of J\f = 4 super Yang-Mills theory on S" 1 x S 3 , when 
restricted to commuting matrices for the homogeneous modes of A and $j. Our 
results cover the whole range of temperatures and couplings for which the one loop 
potential is valid: TR C A -1 . We use R to denote the radius of the spatial S 3 . As 
well as finding the minimal action distribution, which preserves SO (6) R symmetry, 
in Appendix B we exhibit other interesting saddle points of the effective action with 
various patterns of R symmetry breaking. 

In the supersymmetric sector, the logarithmic repulsion between — and subsequent 
condensation of — scalar eigenvalues is crucial for the emergence of a dual spacetime 
geometry in the large N limit [4, 5]. This effect needs to be considered also at finite 
temperature. A central result of our work is that at high temperatures TR = A -1 / 2 , 
the condensate of scalar eigenvalues backreacts sufficiently onto the distribution of 
the A eigenvalues to cause a new phase transition, in the commuting matrix saddle. 
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1.1 Summary of results on the eigenvalue distribution 

At low temperatures we find analytically that the eigenvalues are uniformly dis- 
tributed as S 1 x S 5 . Here the S 1 is in the A direction whereas the S 5 is in the 
$j directions and has a radius that scales like X/R. 

At the deconfinement temperature, TR C = — 1/ log(7 — 4\/3) w 0.38 [7], the 
distribution on S 1 develops a gap. We show analytically that this implies that the 
full distribution acquires topology S 6 . 

At high temperatures, 1 <C TR <C A -1 , we find two candidate saddle points. There 
is an S 5 distribution that can be described analytically and an ellipsoidal distribution 
with topology S 6 that we obtain analytically for TR <C A -1 / 2 and numerically at 
higher temperatures. The S 5 solution is at a point in the A direction. Both solutions 
preserve the full SO(6) R symmetry. 

At weak coupling, we find that the ellipsoidal solution has lowest action at tem- 
peratures TR C < TR < A -1 / 2 . The first result of this paper is therefore that, in 
the sector of commuting spatially homogeneous fields, the low and high temperature 
phases of weakly coupled J\f = 4 SYM theory on S 3 are characterised by eigenvalue 
distributions with differing topology: S 1 x S 5 versus S 6 . As the temperature is 
increased to the value TR = A -1 / 2 at fixed weak coupling, or equivalently, as the 
coupling is increased up to A = 1/(TR) 2 at fixed high temperatures, we find com- 
pelling evidence for a second order phase transition. In this transition the S 6 solution 
smoothly collapses to S 5 . We summarize the behaviour of these saddles in figure 1 
below. Details will be discussed below. 

2 The effective potential 

Our starting point is the one loop effective potential of J\f = 4 SYM theory on S 1 xS 3 , 
restricted to a certain sector of the theory. We denote the radius of the S 3 by R and 
P — 1/T will be the circumference of the thermal S 1 . The effective potential is the 
result of integrating out all inhomogeneous modes on the three sphere and depends 
only on the homogeneous modes. The latter are the vacuum expectation values of 
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Figure 1: Topology of the dominant eigenvalue distribution in which the only the homo- 
geneous modes of {Aq, <£j} condense and furthermore commute. 

the six scalar fields in the theory, $j, J = 1, 2, . . . , 6, 1 and the time component of 
the gauge field, A . This potential was computed in [18] and also discussed recently 



An important feature of the effective potential is that it is evaluated on the space of 
mutually commuting scalar homogeneous modes. The potential is therefore a function 
of the eigenvalues of the adjoint scalars, which we denote by {<fi p }, p = 1, 2, . . . , N, 
and the eigenvalues of /3A , which we denote by {0 P }. 

It is natural to consider the effective potential on the space of commuting scalar 
field expectation values. In flat space and at zero temperature, these characterize the 
moduli space of vacua of the M = 4 theory, wherein off diagonal modes obtain masses 
\4> pq \- 2 There is no moduli space on S s due to the classical conformal coupling to 
the background scalar curvature. However, around a homogeneous background with 
commuting VEVs for the scalars, off diagonal fluctuations of the scalar zero modes 



obtain a positive mass given by y R~ 2 + \<fi pq \ 2 and those of A have masses \<f) pq \. 
For this reason all the configurations that we consider will be locally stable against 
small fluctuations involving the off-diagonal entries of the scalar field matrices. This 

1 We denote a six vector with a boldface symbol, e.g. <I> and 4>. 



in [19]. 




2 We set |0| = \fY,j4> 2 j and 4> pq =4> v ~ <t> q - 
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is also manifest in the reality of the one loop effective potential below, on the space 
of simultaneously diagonal configurations. 

We should emphasize that whereas the truncation to commuting matrices is con- 
sistent, it will not describe the absolute minimum of the action [17]. We will see 
however, that these commuting saddles display some interesting dynamics. Through- 
out this work, we are only considering saddles in which the matrices commute. 

A second consistent truncation we have made is to consider condensates for only 
the homogeneous modes. At infinite N, the inhomogeneous modes can also condense 
without breaking the spatial 5*0(4) rotational invariance. The full dynamics should 
include these modes also. 

There are three contributions to the potential: the classical conformal mass term 
for the scalars and then two determinants arising from integrating out the bosons and 
fermions in the theory 



s eS [d> p ,e p ] = s^ + s^ + s^ 



The classical term is 



N 



A 



(2-1) 



(2.2) 



P =i 



Integrating out the gauge fields, scalar and ghost fluctuations yields the bosonic 
contribution 
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The fermion fluctuation determinants contribute 
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(2.5) 



The validity of the one loop potential requires weak 't Hooft coupling, A 1, but 
also A <C \/{TR). This is because at TR ~ A -1 the spatial S 3 is the same size as 
the nonperturbatively generated magnetic screening scale. Our objective is now to 
minimise this potential and analyse how the ground state depends on the temperature 
and coupling. The Casimir terms in (2.3) and (2.5) 3 can consistently be ignored, 
since expanding them in gives terms of the form f3R\<p\ 2 (l + 0(R 2 4> 2 )) which, at 
weak coupling, is always subleading compared with the tree level mass term. It is 
important, however, that we do not ignore the <j> dependence of the exponential 
terms because these are enhanced at high temperature TR — > oo, and will play a 
significant role in our analysis in this regime. 



1 t 

3 Low temperature distribution: S x 



The action at low RT < 1 is 



S T R«1 = 7 Yl I^pI 2 ~ Yl bg 



^ p i IX/ I 



sinh 



-N 2 \og2. (3.1) 



This expression neglects terms in the action (2.1) that are exponentially suppressed at 
low temperatures. The logarithmic repulsive force generated between the eigenvalues 
is a generalised Vandermonde type interaction resulting from integrating out the off- 
diagonal fluctuations; its origin can be traced to the ghost and gauge fixing terms 
in the action [18]. At A = 0, when the scalars identically vanish, it reduces to the 
measure factor in the unitary matrix model of [7]. 

The pairwise eigenvalue repulsion is only countered by an external quadratic at- 
tractive potential along the scalar field directions. We should therefore expect the 
lowest energy state to be localized in the scalar directions and maximally spread out, 
i.e. uniformly distributed, along the S 1 . In this case, using the fact that 



log 



-log2 + q^|-^-e-"^cos(^) , (3.2) 



n 

n=l 



and noticing that the periodic pieces will average to zero for a uniform distribution 



3 These are the first terms in the sum over £, involving square roots. 
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of 9 p s, we have an effective action 



s TR «i = — ^ — 2^ l^pl ~ 2 2^ I^mI • ( 3 - 3 ) 

p=l pg=l 

In fact we will find in the next section that this action is valid all the way up to the 
confinement /deconfinement transition. One way to see this is to note that, as we do 
in the next section, that below the confinement /deconfinfement transition the 8 p s are 
uniformly distributed around S 1 . Hence, the averaging over them has the same effect 
in the action as sending T — > 0. 

It is convenient to introduce the dimensionless variables 



x 



p 



(3<t> p . (3.4) 



At large N we can pass to a continuum limit. The eigenvalues are described by a 
joint distribution in seven dimensions, which satisfies the normalisation condition 

N 

A* 



1 N 

-Y 



J d 6 xd9p(x,9) = 1 , (3.5) 



P =i 

although in the present case p(x, 6) = p{x) only. The equation of motion for x is 

n 2 RT 



A 



-x = n [ d 6 x'p(x') X X ' (3.6) 
Jd \x — X \ 



where the eigenvalue distribution has support in some domain D C R 6 and x E D. 

We can see that (3.6) has a solution for which D is an S 5 of radius r, which is 
easily determined. The density is found to be 

5(\x\-r) 

P( x ) = 27r 4 r5 > ( 3 - 7 ) 

where the radius r is 

^1024 

r " niRT 945 ' ^ 
This solution coincides with the zero temperature S 5 solution of [4] which is what 
we should expect in the low temperature limit where thermal effects are negligible. 
Importantly, the smooth S 5 is a consequence of the large N limit. The SO (6) sym- 
metry of the M = 4 theory is left unbroken in this saddle by a global 5*0(6) rotation 
on each of the x p combined with the action of the Weyl group of SU(N) which acts 
through permutations of the x p . 



Let us determine the action associated to this S 5 x S 1 geometry of eigenvalues 

S s5xSl = 2N2j ^ TR J d e x ^ x 2 _ An 2 N 2 J d 6 x J d&x l p ( a; ) p ( a .') l og \ X - X '\. 

(3.9) 

To this order in perturbation theory, we find that the action for this configuration is 

1 q A 10242 1 mm 

The action is lower than that of both the solution in which the scalars have zero ex- 
pectation value, and also of the band solution found in [19], in which all x p acquired 
non-zero expectation values along the same direction in 5*0(6) space, thus only pre- 
serving an SO (5) subgroup of the global symmetry group. It seems very likely that 
this maximally symmetric saddle is the absolute minimum of the effective potential 
at low temperatures, within the sector of the theory we are considering. We prove 
stability against R symmetry breaking perturbations in Appendix A. Furthermore, in 
Appendix B we find various other saddle points with reduced R symmetry and find 
that they have a higher action. 

It is striking that the eigenvalue distribution localises to a hypersurface and does 
not spread out in all the six noncompact dimensions. In fact, we can exclude six 
dimensional solutions to the effective equation of motion for the scalars (assuming 
a uniform distribution on S 1 ), by repeated application of the V operator in six 
dimensions. In particular, acting on equation (3.6) with the operator V 2 (V-), we 
find J d 6 x' p(x')/\x — x'\ z = which implies a vanishing density in six dimensions. 
This then implies that the eigenvalues at large N must be constrained to lie on a 
hypersurface in six dimensions and this is indeed what we find. An argument along 
similar lines was shown in [4], implying that eigenvalue distributions arising from 
similar commuting matrix models were always singular. Localisation to a hypersurface 
is the simplest way to achieve this. 

Note that our analysis is consistent with [7] when A = 0: the emergent S 5 disap- 
pears and only the configuration space of eigenvalues of the Polyakov loop remains. 
As advocated in [2, 4] the large N vacuum configuration of eigenvalues should be 
interpreted as the emergence, at weak 't Hooft coupling in our case, of the S 5 in the 
dual AdS 5 x S 5 geometry. The connection to a dual geometry is discussed in more 
detail below. 
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4 The deconfinement transition: S 1 x S 5 —> S 



This section describes the distribution for TR <C A -1 / 2 . In particular, we will be in- 
terested in the 'deconfinement' or 'Hagedorn' transition that occurs at TR C pa 0.38 [7]. 
We will see that the eigenvalue distribution undergoes a topology change: S 1 x S 5 — > 
S 6 . We will comment on possible spacetime interpretations of this transition in a 
later section. The details of the transition are accessible within standard perturba- 
tion theory in A which is valid for all temperatures TR <C A -1 / 2 . 

To analyze the vicinity of the phase transition we make use of certain properties of 
the equilibrium distribution of the scalar fields at weak coupling. At an equilibrium, 
the tree level attractive potential balances out the pairwise repulsive forces generated 
at one loop. It is straightforward to see that the one loop contribution to the effective 
action (2.1) can only compete with the tree level term if \x\ ~ 0{\). Hence, in 
order to be consistent with the perturbative expansion we should expand the one 
loop contribution in powers of x. This has the effect of reorganizing the pertubative 
expansion. The only caveat is that when the temperature is very high, TR > A -1 / 2 , 
a further reorganization occurs which we shall describe in a later section. In this 
section we shall assume that TR <C A -1 / 2 . 

We now define x = \x and expand the action up to 0(A) to obtain 

S = \°S(°\0) + A (SW(0, x) + S 2Aoop (9)) + 0{\ 2 ) , (4.1) 

where 

N oo 

= E E - V - z B{e- n/TR ) - (-ir +1 z F (e- n/TR )} cos(n9 pq ) . (4.2) 

77, L 

pq=l n=l 

This is precisely the A = effective potential of [7]. The single particle partition 
functions zb and zp will be given below. At the Hagedorn transition, the coefficient 
of the lowest cosine term becomes negative [7]. The one loop term in the action (4.1), 
using (3.2), is given by: 

N j N oo 

s (1) [e,x] = Nti 2 trJ2\x p \ 2 - g E 1 + 2E cos ( n M • ( 4 - 3 ) 

p=l P<7=1 71=1 

It is clear from the form of the expansion (4.1) that to leading order the 9 p s are 
unaffected by the scalars. In other words the 9 P , and their density in the large N 
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limit, behave to leading order exactly as they did with zero scalar VEVs. The scalar 
VEVs themselves are then, to leading order, determined by finding the minimum 
of S {1 \9,x) with the given values for 9 P . Notice that although this term is of the 
same order as a two loop contribution, the latter contribution only involves the 9 P , 
as indicated, and therefore doesn't contribute to the leading x distribution. 

The term in square brackets in (4.3) is simply the delta function n5(9 pg ), restricted 
to even functions. Note that the 9 eigenvalue density is indeed an even function [7]. 
Taking the large N limit and describing the eigenvalues by the joint density p{9, x), 
as previously, we have 

A_s w = 7r 2 TR J d9d 6 Zp(9,x)\x\ 2 -ix J d9d 6 x d 6 i' p(9,x)p(9,x')\x - x'\ . (4.4) 

We now proceed to minimise this action. 
The fact that the reduced density 

m = j*w,*), (4.5) 

is determined by the zero scalar VEV problem, together with preserving SO(6) R 
symmetry, implies that the joint density is 

e^mtm t , («) 

\x\ 5 y/l + {dr/d6) 2 Vo\S 5 

where the unknown function r{9) determines the size of the S b as it is fibred over 
the support of p(9). We recall that that we have already argued that other 5*0(6) 
symmetric configurations for which the scalar field spectral density has a smooth six 
dimensional support are not allowed, as in [4]. From (4.6) we have 

-L,S (1) = n 2 RT J d9p(9)r(9) 2 - nV2C J d9 p(9) 2 r(9) , (4.7) 

where 

C= « (4. 8) 
945tt v ; 

Completing the square gives 



The final term only depends on 9 and simply contributes to the two loop order 
distribution of the 9 p s. Therefore it can be ignored for our purposes. However, it 
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will potentially influence higher loop computations for determining the order(s) of 
possible phase transitions at finite coupling [7] 
Hence, for a minimum we have 



r[e) = ^ {e) ■ | (4 ' 10) 

This result directly connects the shape of the eigenvalue distribution and the eigen- 
value density, and implies that the topology change S 5 x S 1 — > S 6 occurs as T is 
increased through the phase transition. This is because p{9) changes from the uni- 
form to a gapped distribution with support —9 < 9 < 9 , where p(±6*o) = [7]. 
When p(9) is uniform, the distribution is S 1 x 5 5 , as in the previous section. Once 
p(9) becomes gapped, then the distribution is an S 5 fibred over an interval with the 
size of the S b vanishing at the endpoints. This is topologically an S 6 . 

We can now check the consistency of our expansion in x by investigating whether 
the higher terms in the expansion are finite on the solution (4.10). The potentially 
dangerous higher terms, coming from expanding e~ n ^ Xp ^ cos(n9 pq ), have the form 

N oo 

1 nx pq | m cos(n9 pq ) 

pq=l n=l 

~ f di sin 4 i f° d9 \p(9)p(9') (p(9) 2 + p(9'f - 2p(9)p(9') cos (m+1)/2 " 

JO J -9 1 J 

(4.11) 

Here, we used the fact that 

oo 

n m cos(n9) = nd r g n 5(9) , (4.12) 

n=l 

and integrated by parts. The question is whether the terms (4.11) are finite. The 
potential problem occurs at the edges of the distribution 9 = ±9 where there are 
possible singularities. We know [7], see also the following section, that in the vicinity 
of the edge 9 = 9 the density behaves as p(9) ~ ^9 — 9. It is not difficult to show 
that the integrals are completely regular at 9 = 9 . Hence the analysis is consistent 
across the transition and into the high temperature phase (assuming (TR) 2 <ti 1/A). 

In summary, inclusion of the lowest order corrections in A leads to an interpretation 
of the Hagedorn/deconfinement transition as a topology changing transition, S 1 x 
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S 5 — > S 6 , in eigenvalue space. The possible implications of this geometric transition 
for the dual spacetime will be discussed below. 

5 Intermediate temperatures: S ellipsoid 

We now determine the joint eigenvalue distribution in the temperature range 1 <C 
TR <C A" 1 / 2 and provide further evidence for the appearance of the S 6 topology 
above the Hagedorn/deconfinement temperature. In this region, the 9 P have the 
Wigner semi-circular distribution. To see this we have to analyse for TR ^> 1. 
First of all, recall that the single particle partition functions in (4.2) are given by [7] 

, , 2x(3 + Qx - x 2 ) . . lQx 3 / 2 

*b(x) = {1 _ x)3 , z F (x) = ^—^ . (5.1) 

For TR > 1 this implies 

z B {x) -> 16(TR) 3 /n 3 , z F (x) -> lQ(TR) 3 /n 3 , (5.2) 

and we expect that the 9 P are small at high temperatures, in which case using 
X^°=i(2™ — l) -2 = 7r 2 /8, we have in this limit 

N N N 

<?(°) - - (log|^|-27r 2 (Ti?) 3 ^) = - ^log|^|+4iVvr 2 (Ti?) 3 ^^ 2 . (5.3) 

pq=l P9=l p=l 

This is the action of the conventional Hermitian matrix model with quadratic poten- 
tial. The saddle point equation is 

N 1 

E T ^ = ±Nk\TR)% (5.4) 

q&p)=l P q 

In the large N limit we can solve for the density in the usual way, by introducing the 
resolvent 

1 N 1 

«w = j?E^. < 5 - 5 > 

in terms of which (5.4) becomes 

u(6 + ie) +uj(6-ie) = 8tt 2 (TR) 3 6 , (5.6) 
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for — 6q < 9 < 6q and e is an infinitesimal which imposes on the left hand side a 
principal value. It follows from (5.5) that 

uj(6 + ie) - uj(6 - ie) = -2mp(9) , (5.7) 

where again — 9 < 9 < 9 . It then follows from (5.6) and (5.7) that the resolvent 
uj(x) is an analytic function of x which has a square root branch cut between x = ±9o. 
Since u(x) must go to zero for large |x| this determines uniquely 



u(x) = An 2 {TRf [x - ^Jx 2 - . (5.8) 
By taking the discontinuity across the cut we obtain 



p(9)=Av{TRf^9l-9*. (5.9) 
Finally, normalising the density gives 

ol = ^fW ■ (5 ' 10) 

Using (4.10) we see that in this region of intermediately high temperatures the 
combined density has support on an ellipsoid with topology S 6 given by 



-x 2 ■ ■2;r [ TI!f9- 1 . 



\UC 2 TR 



(5.11) 



Given that TR\^X <^ 1, this ellipsoid is very elongated in the 9 direction. The 
eigenvalue density (4.6) is largest at the highly pointed tips at x — 0. 



6 High temperature distributions: S Q versus S 5 

This section studies the distribution at high temperatures 1 TR <^ 1/A, which 
includes TR ~ A^ 1 / 2 . We will see that in this regime the VEVs of the scalar fields 
are no longer determined by the 9 P , but rather the full coupled system must be 
considered. The new phenomenon we will find in this regime therefore could not be 
seen in previous analysis that neglected these scalar VEVs. 

It is worth remarking that for TR > A -1 / 2 , the S 3 size exceeds the electric or 
Debye scale, and perturbation theory in A is replaced by a perturbation theory in \f\ 
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due to infrared effects. Importantly for us, our one loop effective potential remains 
unaltered by these resummations since the associated momentum integrals at high 
temperatures happen to be insensitive to IR effects. 

The analysis of the previous section, however, will break down at high tempera- 
tures where TR > A -1 / 2 for a different reason. At these temperatures, the distribution 
of the x will begin to affect the 9s. This is precisely the region when the S 6 begins 
to look spherical in the dimensionless variables x and 9. However it is important to 
realise that the analysis in this section has an overlap with the analysis of the last 
section when 1 < TR < A~ 1/2 . 

To see what goes wrong with the earlier analysis look at the corrections (4.11). 
These are 

N oo ^ N 

^ -\nx pq \ m+1 cos(n9 pq ) = vr ^ \x pq \ m+1 d™5(9 pq ) . (6.1) 

pq=l n=l P9=l 

We can estimate the behaviour of these terms as a function of A and TR by using the 
fact that from (5.11) when TR >lwe have x ~ \(TRf/ 2 and 9 ~ (TR)' 3 / 2 . Hence, 
the correction goes like (A(T R) 2 ) m+1 . Clearly these corrections cannot be ignored in 
the high temperature regime. 

However, we expect that when TR ^> 1 then the eigenvalue distribution will 
satisfy 9 p ,\x p \ <C 1. With these assumptions, the high temperature action is found 
to be 

Stb>i = £W 2 4E log (|« M | 2 + <&) + - W £ (|^| 2 + 29 2 pq ) . 

p=l P9=l P<7=1 

(6.2) 

The validity of the assumptions is verified a posteriori from the solution. The large 
iV equations of motion following from this action may be written 



Px= [ (fx' d9' p(x', 9') ^ 

J \x — x' 



x — X 



2 + (9 - 9') 2 

Q9 = j d«x> d9> P (x>, 9') ][X _Jf^[ e _ dr , (6-3) 

where P = 7i 2 TR(l/\ + 2R 2 T 2 ) and Q = 4tt 2 R 3 T 3 . These equations have more than 
one interesting solution. Once again, by application of the seven dimensional operator 
V 2 (V-) in the (x,9) space, we can argue that non-trivial solutions to the equations 
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of motion must be hypersurfaces in seven dimensions. For a given P and Q there 
are two solutions which preserve the SO(6)r symmetry. One is topologically S 6 and 
the other is topologically S 5 . We expect one of these maximally symmetric solutions 
to have the lowest action. Before discussing the solutions, we comment on a scaling 
property of the solutions with P and Q. In appendix C we present various solutions 
that do not preserve the full R symmetry. 

6.1 Scaling of the solutions and action 

The effect of changing the coupling or the temperature of the theory is incorporated 
in the values of P and Q in the high temperature equations of motion (6.3). Under 
the rescaling 

P = fiP, Q = fiQ, (6.4) 

the equations of motion and the normalisation condition (3.5) imply that the solution 
remains the same and scales to 

x = — x, e = —e. (6.5) 

The action evaluated on the scaled solution becomes 

S[x,9]=S[x,e] + ±lo gf i. (6.6) 

We can use this scaling property to fix Q and work with P/Q. When we do 
numerics shortly, it will be convenient to work with Q — 1. Even though this value is 
not in the high temperature regime, we simply have to scale the solutions as we have 
just described. Note that 

P = 1 1 1 , fi 7 ^ 

Q A{TRf\ 2' 1 ' 

so that P/Q — > oo corresponds to (TR) 2 X <^ 1, whereas P/Q — > 1/2 corresponds to 

high temperatures (TR) 2 \ 3> 1 with the caveat (TR)\ <^ 1. The solution in this 

regime only depends on the combination (TR) 2 X and so, increasing or decreasing this 

parameter can also be interpreted as varying the value of the 't Hooft coupling A, at 

a fixed temperature and at weak coupling. 
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6.2 S solutions 



These have 9 = and preserve the full R symmetry group of the theory in the large 
N limit. We can write the eigenvalue density 

where the radius is 

' (6-9) 



'IP 

Evaluating the action on the solution gives 

W2 S ss = Y4 + l l ° g2P - (6 ' 10) 
These solutions are fully collapsed in the 9 direction, p(9) = 5(9). At low tempera- 
tures, these are unstable saddle point configurations, but play an important role at 
high temperatures as we show below. 

6.3 Ellipsoidal S 6 solutions 

We now find the solutions which preserve SO(6)r symmetry and where the eigen- 
values spread out into a (closed) six dimensional surface in R 7 . These are the high 
temperature continuations of the ellipsoids we found at intermediate temperatures. 
We can parametrise the surface by 

= f(r), (6.11) 

where as before (although now without the factor or A) r = \x\. The induced metric 
on this surface is 

ds 2 = (1 + f{rf)dr 2 + r 2 (# 2 + sin 2 dQ 2 4 ). (6.12) 



Thus, for instance, for a sphere of radius A we have f(r) = ±V A 2 — r 2 while for an 



ellipsoid with semi- axes A and B, f(r) = ±B^l — r 2 /A 2 . Note that the parametri- 
sation we are using forces us to consider the solution in terms of two separate com- 
ponents. We now proceed to derive an action whose variation determines f(r). 

By 5*0(6) symmetry, the most general spectral density for the eigenvalue distri- 
bution at large iV must be of the form 

p(x, 9) = 9 (r) (5(9 - f(r)) + 5(9 + f(r)) (6.13) 
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where f(r) is taken to be positive. We can now determine the action restricted to 
the SO (6) symmetric ansatz. As we will be reducing the problem to one dimension, 
it is convenient to introduce the effective spectral density 



G(r) = 27r 3 rVi + f'{r) 2 g{r) . (6.14) 



The high temperature action is therefore rewritten as 



J drG(r)(Pr 2 + Qf(r) 2 ) - £- j dr j dr' d0sin 4 (6.15) 

x G(r) G(r>) log (|r - rf + (f(r) - fir')) 2 ) (|r - rf + (f(r) + f{r')f) . 

We should also add a Lagrange multiplier term that implements the normalisation 
constraint: fi (J drG(r) — l). By functionally differentiating with respect to G and / 
we can show that this action consistently leads to the equations of motion restricted 
to our 5*0(6) ansatz. 

It is possible to perform the <fi integral in (6.15) using contour integration. The 
answer is 



7 

log rr' + — - log 2 



j^S = J drG(r)(Pr 2 + Qf(r) 2 )-^ J dr J dr'G(r)G(r') 

+ log (v^TT + VK^T) - ^ 3^v^l + 4(^^-1) 1 

+ n [J drG(r) - lj , (6.16) 



where 

R 2 = r 2 + r' 2 + (f(r)-f(r')) 2 

2rr' ' 1 ' 1 

In (6.16) we have reduced the problem to one dimension. Variation of this action 

leads to a pair of integral equations for f(r) and G(r). We have solved these equations 

numerically using a simple Monte Carlo algorithm. The algorithm discretises the r 

axis into iV points: f drG{r) — > Y^i=i an d then minimises the multiparticle action 

(6.16) by relaxation. A test of our code is that at the point with enhanced SO (7) 

symmetry, when P — Q, it correctly reproduces a six sphere of radius 1/V2P. The 

action of this configuration may be found analytically to be 

1 



N 2 



S S 6 



107 1, P * i 

+ - log-. (6.18) 



p=q 120 2 °2 
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The numerically obtained action agrees with this exact result to three significant 
figures, even with a fairly small number of points, N ~ 60 or so. Our main interest 
is for P ^ Q, for which we currently do not have analytic results. 

The equations that we are solving numerically are very similar to those consid- 
ered in the recent paper [5], which was in a supersymmetric context. The two main 
differences are firstly that we are looking for the semiclassical, N — > oo, saddle point, 
whereas [5] perform a quantum Monte Carlo simulation to study the finite N wave- 
function. Secondly, we have restricted to an 5*0(6) invariant ansatz before applying 
numerics. This gives a more efficient use of eigenvalues, which are diluted in fewer 
dimensions. However, it means that the SO(6) symmetry is a (consistent) input as- 
sumption. The numerics in [5] did not assume this symmetry but found it in their 
results. This provides the complementary information that the 5*0(6) ansatz is stable 
against R symmetry breaking perturbations. 

Figure 2 shows the eigenvalue distributions obtained for three values of P/Q. The 
form is roughly what we could have anticipated, given that as P is increased the 
external force pushing the eigenvalues to small r becomes greater. The eigenvalue 
distributions appear to be ellipsoids to a very high degree of accuracy. This is per- 
haps not surprising given our results from the intermediate temperature regime. The 
numerics are not accurate enough to reliably read off a density profile. 

The gap in the data points at small r in figure 2 is due to the fact that the effective 
density, G(r), is very low in this region. This is due to the r 5 term in (6.14). Note 
however that there is a data point at r = for all three curves. The solid lines in 
figure 2 are the ellipses determined by the extent of the distribution in the r and 
/ directions. These curves match all the numerically obtained data points to two 
significant figures. 

Given the appearance of ellipsoids, a rescaling-invariant question we can ask is 
how the axis ratio of the ellipsoid varies with P/Q. Figure 3 is a plot showing the 
axis ratio at several values of P/Q. There appears to be a linear relation between 
these quantities to good approximation over the range plotted. Figure 3 also shows 
a linear fit to the data points for P/Q > 2. The best fit line is 

^ = -0.81 + 1.51^, (6.19) 
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Figure 2: Output from the Monte Carlo simulation for {P, Q} = {5, 1}, {2, 1}, {1, 1} with 
N = 500, 500, 1000 points, respectively. The solid lines are ellipses with axis ratio deter- 
mined by the data. 

where a and b are the lengths of the axis in the / and r directions, respectively. 

In the limit where X(TR) 2 <C 1, or P/Q large, we can compare the axis ratio with 
that of the ellipsoid we found in the intermediate temperature regime (5.11). In terms 
of P/Q it follows from (5.11) that the axis ratio for those ellipsoids is also linear in 



This is fairly close to (6.19), especially given that the fit (6.19) included a range of 
points outside the region of validity of (5.11). A more precise matching comes from 
taking, for instance, the value P/Q = 20. This corresponds to X(TR) 2 ~ 1/80, and 
so should be well described by the intermediate temperature analysis. For this value 
we find that the numerical axis ratio is 29.6 whereas the analytic result is 29.0. This 
close agreement is a test of both the numerics and analytic results. 

Finally, we can compute the action of the solution as a function of P/Q. The 
result is shown in figure 4 for Q — 1. We have also included in the plot the action of 
the round S 5 solution and the action of the two dimensional ellipse solution discussed 



As with the axis ratio, at large P/Q we can compare the numerically obtained 



P/Q 



a _ 945vr P 
b ~ 2048 Q 




(6.20) 



in [19]. 
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Figure 3: The axis ratio of the ellipsoidal distribution as a function of P/Q, with Q = 1. 
The data points are from the Monte Carlo numerics whereas the solid line is a linear fit. 

actions with the analytic result in the intermediate temperature regime. At P/Q = 20 
the numerically found action is 1.04 and the action of the ellipsoid (5.11) is 1.09. 

We see that throughout the range plotted, P/Q > 1, the ellipsoidal eigenvalue 
distribution has the lowest action and therefore dominates the partition function. 
Therefore, at sufficiently weak coupling, the ellipsoidal distributions describe the vac- 
uum of the theory. However, as P/Q is lowered past 1, recall that lowering P/Q 
corresponds to either increasing the temperature at fixed weak coupling or increasing 
the (weak) coupling at fixed temperature, it looks like the curves for S 5 and the six 
dimensional ellipsoid might cross. This raises the prospect of a new phase transition 
in the high temperature regime as a function of coupling. 

7 A second order phase transition: S — > S 

In figure 4 we saw that as the coupling is increased at fixed high temperature, the ac- 
tion for the S b solution comes closer and closer to that of the, increasingly anisotropic, 
ellipsoidal S 6 solution. At the point with enhanced SO (7) where we have analytic 
results, P = Q = 1, the difference between the actions is surprisingly small, given 
that P/Q = 1 is not particularly close to the asymptotic regime P/Q — > 1/2 

^ [S* -S sB ] p=q=1 k -0.0098. (7.1) 
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Figure 4: The action of three different eigenvalue configurations as a function of P/Q with 
Q = 1. The data points are from the S 6 ellipsoid. The solid curve is the round S 5 solution. 
The dotted curve is the two dimensional ellipse solution. 

What happens as P/Q is decreased below this value? We now give four arguments 
that, taken together, strongly suggest that a second order phase transition occurs in 
the commuting saddle at X(TR) 2 = 1, or P/Q = 3/4. This is the curve that we 
included in figure 1 above. 

7.1 The actions converge 

Figure 5 shows the action from the numerical computation together with the ana- 
lytically computed action of the round S 5 solution (6.10) as a function of P/Q, for 
0.7 < P/Q < 1.2. We see that to within the accuracy of the numerical computation, 
the error seems to be under 0.5%, the numerical data joins the curve of the S 5 action 
somewhere around P/Q ~ 0.75 or 0.8 and then follows this curve for lower P/Q. 

The simplest interpretation of figure 5 is that beyond some critical coupling, such 
that P/Q ~ 0.75 or 0.8, the S 5 saddle becomes the lowest action configuration. The 
Monte Carlo simulation therefore converges on the S 5 saddle below this value. If this 
indeed occurs, it indicates that there is a quantum phase transition at the critical 
coupling in which the vacuum manifold of eigenvalues changes topology: S 6 — > S 5 . 
It is a quantum phase transition in the sense that it is not associated with any 
symmetry breaking and the different phases have topological characterisations. It 
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Figure 5: The data points are the numerically computed action (with N = 150). The 
curve is the action of the S 5 solution. 

is possible however to write down non-topological "order parameters" that go from 
non-zero to zero as we cross the transition, such as (x 2 ) 2 — (x 4 ). It should be clear 
from this order parameter that the transition will not be visible if the scalars are 
integrated out. 

The convergence of actions does not tell us whether the possible transition is first 
or second order. We will now show analytically that at P/Q = 0.75 the S 6 saddle 
becomes completely localized in the 9 direction and becomes a round S 5 . Further- 
more, we we will find zero modes about the S 5 solution at precisely P/Q = § which 
provide compelling evidence that there is indeed a second order phase transition. 

7.2 The S 6 collapses to S 5 

Motivated by the numerical results above, we now look for an analytical condition 
on the parameters P and Q, for the S 6 solution to collapse to an S 5 . The central 
fact we will use is the existence of the S 6 saddle for an appropriate parameter range, 
a fact which we have established both numerically and analytically for a wide range 
of temperatures. In the temperature range TR 3> 1, the S 6 topology solves the 
equations (6.3). 

Since the S 6 solution is actually a round S 5 fibred over a finite interval in the 9 
direction, we may differentiate the second equation of motion (6.3) with respect to 9 
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to obtain a condition which must be fulfilled by the same solution : 

Q = f (fx' dO' p(x', 9') \x-x'\ 2 -{0-6'f 

J \\x - x'\ 2 + (9 - 9'Y) 2 V ; 

We emphasize that this will be automatically satisfied by the 5*0(6) symmetric S 6 
topology, since the equation of motion holds at each point on the support of p(x, 9) 
in the 9 direction. 

We now want to find at what temperature the extent of the effective distribution 
in 9 shrinks to zero size, and thus the S 6 collapses to an S 5 . In such a limit, we 
expect the density function p(x, 9) to smoothly approach (6.8) with the S 5 radius 
equal to 1/V2P. Substituting these into (7.2), we find 

hi- 

It is easy to check that the integral in (7.2) is well behaved in the limit that 9, 9' 
approach zero, and that the limiting value of | is approached from below. Hence 
the S 6 topology ceases to exist for Q/P > | which translates to TR > A -1 / 2 , and 
at these temperatures with a fixed weak coupling, the equations of motion are only 
solved by the S 5 configuration of eigenvalues. 

We remark that this transition from S 6 to S 5 cannot be seen in the A = theory, 
since the critical temperature is T = l/(\f\R) , and is driven by the presence of the 
scalar expectation values in the J\f = 4 theory. 

7.3 The S saddle develops a zero mode 

Now consider the stability of the S 5 solution (6.9). As mentioned above, we expect the 
S 5 solution will be stable against perturbations in the x p directions that would break 
R symmetry. It is important to check this and we do so in Appendix B. However, 
something interesting occurs if we consider perturbations in 9 P . 

Returning for the moment to the discrete system, the quadratic action for fluctu- 
ations 89 p is given by 



p=l P<?=1 



P'l\ 



4 We may do this consistently for every value of 9 lying inside the distribution, and by continuity 
we can aslo apply the condition to all points approaching the edges of the distribution in 0-space. 
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It is straightforwardly shown that in the large N limit, the N eigenvalues are all equal 
to N(Q — |P). Hence, the solution is only stable when Q/P > |, that is 

TR > -L . (7.5) 

When Q/P < | it is unstable to spreading out in the 9 direction. We have already 
found the high temperature configuration with 9 spread out: it is the S 6 ellipsoid. 

The appearance of zero modes in the S 5 vacuum and negative modes for TR < 
A -1 / 2 , in conjunction with our earlier observation that the S 6 solution merges with 
the S 5 at TR = A -1 / 2 , provides clear and solid evidence for a second order phase 
transition at this temperature. 

We remark that although there are N zero modes at the transition there is only one 
S 6 vacuum due to the action of the Weyl group permutations on the joint eigenvalues 

(4> P ,e P ). 

7.4 The eigenvalue distributions converge 

In further support of this picture, let us look in more detail at the output of the 
numerics as P/Q — > f. The output no longer resembles the ellipses shown in figure 
1 but rather gives a cluster of points around a radius that is very close to the radius 
of the S 5 , l/v / 2~P as found in (6.9) above, and / = 0. This is what we would expect 
to see if the density of eigenvalues of the S 6 were accumulating on the equator as 
the ellipsoid becomes increasingly squashed. In the limit of this process, the solution 
becomes an S 5 . 

The three graphs in figure 6 illustrate this process. As P/Q — > § , the eigenvalues 



start to cluster around 0.81 ~ 1/V2 x 0.75. As in the previous section, the gaps in 
the distribution are partly due to a low effective eigenvalue density (6.14) as r — > 0. 
However, we clearly see the gap growing as the eigenvalues become more densely 
clustered around the equator of the squashed ellipsoid. 



8 Dual spacetime interpretation 

A striking recent development in the supersymmetric sector of the zero temperature 
AdS/CFT correspondence is that the scalar field eigenvalue distribution can be di- 
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Figure 6: From left to right, the eigenvalue distributions with P/Q = 0.9,0.8 and 0.75. 
Numerically computed with N = 150 eigenvalues and Q = 1. 

rectly connected to the dual spacetime geometry [2, 3, 4, 5]. The most immediate 
manifestation of this connection is that the S 5 part of the eigenvalue distribution is 
to be identified with the S b part of the dual AdS 5 x S 5 geometry [4]. One can think 
of the eigenvalue distribution as signalling the fact that the spacetime geometry has 
undergone a geometric transition, in which a noncontractible S 5 appears, due to the 
gravitational backreaction of the N D3 branes. More generally, the eigenvalue distri- 
bution is to be identified with the locus in the dual spacetime where the S 3 of the 
conformal boundary geometry degenerates [4, 5]. In AdS?, this is simply the origin 
of the AdS space. 

The question is then whether, in our finite temperature, non supersymmetric and 
weak coupling setup, we can reinterpret the various topology changes undergone by 
our eigenvalue distributions as topology changes in a dual spacetime, analogous to 
the Hawking-Page transition. Although our saddles are not generically the dominant 
saddles, the assumption of commuting matrices does mean that they are the most 
'geometric' saddles at weak coupling. At low temperatures, our S 1 x S 5 distribution 
looks promising. This is indeed the topology of Euclidean thermal AdS at the origin 
of the AdS space. 

The S 1 appearing in the low temperature S 1 x S 5 eigenvalue distribution should 
in fact be associated with the T-dual of the thermal circle in spacetime, insofar as 
T-dualising along a thermal circle is well defined. This is because eigenvalues of the 
Wilson line around a circle, in our case the Polyakov- Wilson loop around the thermal 
circle, yield the locations of D branes in the T-dual transverse circle. Supporting this 
observation, the radius of the S 1 circle is proportional to the temperature T, if we 



26 



renormalize the low temperature S 5 radius (3.8) to unity. 

The S 6 eigenvalue distribution appearing above the first order Hagedorn transi- 
tion has a natural explanation if one takes the notion of T-duality along the thermal 
circle seriously. In this picture, the Polyakov loop eigenvalues represent positions of 
D2 branes, T-dual to the original D3 branes, along the T-dual circle. The gapped 
distribution of the Polyakov loop eigenvalues above the Hagedorn transition is then a 
localized distribution of D2 branes on this circle. A similar localization of D2 branes 
on a spatial circle, at finite temperature, has been argued [20] to produce a near 
horizon geometry containing a non-contactible S 6 . In that case, somewhat tantaliz- 
ingly, an S 1 x S 5 — > S 6 topological transition of the Gregory-Laflamme type was 
predicted from supergravity, and the localized D2 brane configuration was reinter- 
preted as winding mode condensation on the original circle. The interpretation of a 
nonuniform eigenvalue distribution as smeared branes has also been used previously 
to good effect in studying a field theory dual of the Gragory-Laflamme black string 
instability in [21]. 

From the weakly coupled field theory perspective, the connection to a D2 brane 
setup above the Hagedorn/deconfinement temperature and below T = l/(y/\R), 
perhaps suggests that the field theory is in an effectively three dimensional phase 
in that region of the phase diagram. The potential appearance of an S 6 topology 
in the dual spacetime would be consistent with this line of thinking: At any fixed 
temperature, as the 't Hooft coupling is decreased so that string corrections become 
substantial, at a critical coupling, when the size of the thermal circle in string units 
becomes sufficiently small, the theory undergoes a stringy transition resulting in the 
S 6 topology. 

More interesting is the possible dual spacetime interpretation of the S 5 topology 
above the second order transition at TR = A -1 / 2 . It is tempting to identify this 
topology as the deep interior of the big black hole in AdS^ x S 5 wherein the thermal 
circle has shrunk to zero size at the horizon. The simplest phase diagram implied by 
the new phase transition, shown in figure 1, also suggests that this phase might be 
the continuation to weak coupling of the big black hole in AdS space. 

It is natural to try to associate the disappearance of the product S l factor in 
the eigenvalue distribution with the appearance of a black hole in the dual geometry. 
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However, it is difficult to make this connection precise. For a start, in black hole 
geometries the spatial S 3 never collapses. The Euclidean thermal circle degenerates 
at the horizon instead. The zero temperature connection of the eigenvalue distribution 
to the locus of vanishing S 3 s in the dual spacetime will need to be modified in the 
non-supersymmetric thermal situation for a spacetime interpretation of the eigenvalue 
distributions to be possible. Since these distributions characterize the ground state 
or the deep IR of the thermal field theory on S 3 , it is possible that they provide 
information about the deep interior of the dual geometry. 

9 Discussion and conclusions 

In this paper we have studied weakly coupled finite temperature Af = 4 super Yang- 
Mills theory on a spatial S 3 , in the large N limit. The main difference between 
our work and previous treatments, [7] being the most important, is that we have 
considered the effects of condensed eigenvalues for the six scalar fields of the theory, 
as well as the time component of the gauge field, A . These scalar eigenvalues have 
been studied recently in the zero temperature theory [4, 5] where their condensate 
is directly related to the dual spacetime geometry. The scalar eigenvalues condense 
despite having conformal and thermal masses because there is a logarithmic repulsion 
between the eigenvalues, which at large N overcomes the mass squared terms. 

We have only studied certain special saddle points of the full large N theory, in 
which the homogeneous modes of the fields {A , $j} condense and commute with 
each other. Although these are not generically the dominant saddles, and thus do not 
determine the phase structure of the theory, they have several interesting features. 
Firstly, they are tractable saddles that preserve the full SO(Q) R symmetry of the the- 
ory. We have seen that they undergo nontrivial dynamics as a function of temperature 
and coupling. More speculatively, because these are the most geometric saddles at 
weak coupling, one can wonder whether they have any connection with the geometry 
that arises in the strongly coupled theory via the AdS/CFT duality. For instance, 
if a mechanism similar to that described in [4] is responsible for the emergence of 
spacetime geometry at finite temperature, then the strongly coupled theory will be 
described by a commuting saddle point. In this scenario, our weak coupling saddles 
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would be continuously connected to the strong coupling supergravity geometry. 

At low temperatures and through the Hagedorn phase transition into the de- 
confined phase, the picture of [7] is not fundamentally modified. The eigenvalue 
distribution of A which that paper considered becomes part of a higher dimensional 
distribution: S 5 x S 1 at low temperatures and an ellipsoidal S 6 above the transition. 
However, at sufficiently high temperature TR = A -1 / 2 , or alternatively, when the cou- 
pling becomes sufficiently strong (but still small) we found evidence for a new second 
order phase transition, in the commuting saddle. This transition is directly due to the 
backreaction of the scalar field eigenvalues on the eigenvalues of A . Therefore, this 
transition could not have been seen without inclusion of the scalar field eigenvalues. 
In terms of distributions, at the transition: S 6 — > S 5 . 

We presented four pieces of evidence for a second order phase transition in the 
weakly coupled theory at A = l/(Ti?) 2 , with TR 3> 1. Firstly, we found numerically 
that the actions of the S 5 and S 6 saddles meet at this coupling. Secondly, we showed 
analytically that the S 6 solution collapses to S 5 . Thirdly, the S 5 solution develops 
a zero mode at precisely this point. Finally, we numerically showed how the S 6 
distribution becomes increasingly squashed and approaches an S 5 at this coupling. 

An application of the results in this paper is that they will allow the computation 
of the Maldacena-Polyakov loop at weak coupling in the commuting saddles. This 
is a thermal Wilson loop involving both the gauge potential and the scalar fields, 
which arises naturally in the M = 4 theory. These loops may be computed at strong 
coupling [23] and are therefore an interesting observable to compare at weak and 
strong couplings. It would be interesting to understand how the S 6 — > S 5 phase 
transition we have discussed is reflected in the value of these loops. 
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A Stability of the solutions 

In this appendix, we will consider the stability of some of our solutions. In particular, 
we shall consider both the high and low temperature S 5 solutions. We only consider 
the stability to perturbations in x p and 9 P , and not to other modes that we have 
integrated out. 

The high temperature S 5 

We begin with the high temperature S 5 solution of Section 6.2. By choosing 
^2 P x p = TlipQp = 0' an d scaling the x p and 9 P appropriately, the discrete action (6.2) 
has the form 

N N 

S = ~\ £ log ( + <&) +NJ2(\* P \ 2 + °0l) (A-l) 

pq=l p=l 

with a = Q/P, The S 5 solution corresponds to 9 P = and \x p \ — l/y/2. Expanding 
around the solution to second order in the fluctuations, we have x p = -j^Qp + Sx p , 
where Q p is a unit 6-vector, we have 

u-s± (NMH t (*^-{<b=g^)y ■ ■ . (a,) 

p=l pq=l VI VI 

At large N, let us estimate the orders of of the off-diagonal relative to the diagonal 
terms in the quadratic form of (A. 2). Firstly the off-diagonal terms. One might think 
that when Sfl = Q p — fl q is small these can be large relative the diagonal terms. The 
question is: how small can 5fl be for A^ points distributed uniformly on S 5 ? If it were 
an S 1 then this would be ~ 1/N and so the sum over pairs in (A. 2) would be of order 
A^ 2 and so the off-diagonal terms would dominate and there would be instabilities. 
However, on S 5 the distribution of the relative angle cos i? = Ct p ■ fl q is weighted by 
a factor sin 4 Hence, near t? = the average separation in t? is 8$ where 



— / d$ sin 4 ^ = 1 (A.3) 

OH ' 







30 



i.e., 5$ ~ A r ^ 1//5 . Hence, the off diagonal terms in (A. 2) are at most 0(N 2 ^ 5 ) and 
consequently are subleading at large A" and we can ignore them for the purposes of 
establishing stability. 

To summarize, the issue of stability is determined solely by the diagonal terms; 
however, we shall find it necessary to go up to quartic order in order to settle the 
issue. 

To all orders, we have 

N 

P =i 

N 

- J] log ((i + V2n p ■ 5x p )(i - n p ■ n q ) + \5x p \ 2 + e 2 p - V28 ±Xp • nj) , 

9=1 

(A.4) 

where S±x are the variation perpendicular to fl p , i.e. tangent to S . We can evaluate 
SS by replacing the sum over fl q by an integral f d 5 tt and use 

y^jF(Op-f2 g ,<5j_a;p-r2g) — >■ / dip sin 3 ip / sin 4 i} J r (cos , d,\S±x p \sm'd cosip) . 

q 71 JO JO 

(A.5) 

This gives to quartic order, 
N N 

SS =— y ((3a - 4)6* + 2(O p • 5x p f + 2y/2(Q p ■ 5x p ) (\5x p \ 2 + 20 2 ) 

,\ V (A.6) 



+ \x p f + ^x p \ 2 Q 2 p + §Q^ + 



From this, a careful analysis reveals that the solution is stable so long as a > | as 
claimed in the text. Notice that the fluctuations tangent to the S 5 are only stable to 
quartic order. 

The low temperature S 5 

We can prove the stability of this solution using exactly the same arguments. First 
of all, the action from Section 3 is 



N 



p=l P9=l 

and the solution consists of an S 5 with (f> = rd p , where we have defined the radius 
r = 1024A/(9457T 3 .R). As above only the diagonal terms play an important role in 
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the large N limit. Similar methods give 



5S 



Nir 2 ) 
8R 



\(3 



((Op • W P ? + • ^p)((«p • ^ P ) 2 - 6|^ p | 2 ) 



(A.8) 



p=i 



893025tt 



(5(«p • <%) 4 + 24(O p • ty,)^/ - 24|«S0 p | 4 )) + • • • . 



33554432 



A careful analysis of this reveals expression proves that the solution is stable against 
all fluctuations. 



B A bestiary of saddles breaking R symmetry 



As well as the absolute minima of the action discussed so far, we can also find various 
solutions to the effective action that do not exhibit SO (6) invariance. These solu- 
tions break the maximal R symmetry in certain patterns as we describe below. An 
important observation is that at any temperature they have bigger action than the 
dominant maximally 5(9(6) r symmetric solutions which we have described so far. 
Therefore they are all unstable or metastable. These solutions may be relevant for 
studies of the theory at finite chemical potential. 

One interesting question raised by the existence of these saddle points is whether 
they survive into the strong coupling regime. If they do, then presumably one should 
expect to find corresponding (unstable) supergravity solutions with reduced R sym- 
metry preserved. Of the solutions that we are about to list, the lower dimensional 
spheres were essentially considered in [22], whereas the products and fibrations of 
spheres are new. 

As well as the solutions described here, there are also special two dimensional 
'Coulomb gas' solutions described in [19]. In those solutions only one scalar field is 
non-zero. That solution is also different to those listed below in that it uses the 6 
direction in a non-trivial way. 
Low temperatures: TR <C 1 

Within this range of temperature, the 6 distribution is uniform. Therefore there 
is always an S 1 part of the eigenvalue distribution. In the order of increasing action, 
the solutions we have found are as follows: 
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• S 1 x S 1 x S 3 : The eigenvalue density is 

n(R , 5(\x\ - A)5(\y\ - aA) 

where a is the ratio of the S 3 and the S l and the vectors x and y span 2D and 
4D planes respectively. The equations of motions for x and y determine the 
radius A and the ratio a as 

A=^^, (B.2) 

with C\ = 3.8921 . . . and a = 1.4344 The action evaluated on this solution 

reads 

W2 Ss^xs* - - m 4?r6 • (B.3) 

The difference between the actions of the maximally symmetric S 1 x S 5 solution 
and this S 1 x S 1 x S 3 solution is then 

T F> 1 T Ft 1 

— W 2 AS = — Jp I^x^xs- " ^ xfl «] « 0.98 x 10- 5 . (B.4) 

• S 1 x S 2 x S 2 : The eigenvalue density is p(9,x) = 327r 3^4 ^(l^| — — 
where the vectors x and y span 3D planes. The radius A, action and difference 
in actions with the S 1 x S 5 solution are 

2A(2 v / 2-l) 
15ti 2 RT ' 

J_ A 8 / 2y^2- 1 N " 

TR 1 



A iV 2 



AS « 0.103 x 10~ 4 . (B.5) 



S 1 x S 1 x S 1 x S 1 : The eigenvalue density is p(#, x) = 167r 4^3 ^(l^l ~~ — 
A)<5(|z| —^4) , where the vectors x, y and iTspan 2D planes. The radius A, action 
and difference in actions with the S 1 x S 5 solution are 

A = ~* [C 2 = 98.7075...], 

1 q A3 

-AS « 0.209 x 10" 4 . (B.6) 

A A/ 2 



33 



S 1 x S 4 : The eigenvalue density is p(0,x) = 16 ^ A 4 ^(\^\ ~ A)5(xq), where the 
vectors x span a 5D plane. We have 

12A 
35ir 2 TR ' 
1 „ A 12 2 



Af^ slxS4 i?T35 2 7r 2 ' 

^AS « 0.144 xl0~ 3 . (B.7) 
A 

S 1 x S 1 x S 2 : The eigenvalue density is p(6,x) = 167r3 ^ 2A3 £(|x| — — 
aA)5(a; 6 ) , where a is the ratio of the S 2 and the S 1 and the vectors x and 
span 2D and 3D planes respectively. We have 

A = -^P- [C 3 = 5.4034...], 

1 A (l + a 2 )C! 
— — b oi v ei v s;2 = — — — a = 1.2404... , 

^2 S xS xS RT g47r 4 L J > 

-AS' « 0.165 x 10- 3 . (B.8) 

A A/ 2 

S 1 x S 3 : The eigenvalue density is p(0,x) = 4w l A a 5(\x\ — A)5 2 (y) , where the 
vectors x and y span 4D and 2D planes. We have 

16A 



A = 



15tt 3 RT ' 



1 q A 16 2 



JV 2 SxS ~ RT 15 2 tt 4 ' 

——AS w 0.374 xl0~ 3 . (B.9) 
A I\ 

S 1 x S 1 x 5 1 : The eigenvalue density is p(6, x) = ^s5(\x\ -A)8(\y\ - A)5 2 (z) , 
where the vectors x, y and z span 2D planes. We have 

A = [C 4 = 18.912...], 

1 AC 2 
F 5slxSlxSl = "iff 32^' 

t n 1 

——AS « 0.428 xl0~ 3 . (B.10) 
S 1 x S 2 : The eigenvalue density is p(9,x) = 8 J, A2 8(\x\ — A)5 3 (y) , where the 
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vectors x and y span 3D planes. We have 

A 



3n 2 RT ' 



JV 2VxS2 " Ti^vr 2 ' 

^-^AS ~ 0.796 x 1(T 3 . 

A iV 2 

(B.H) 

S 1 x S 1 : The eigenvalue density is p(9,x) = j^jS(\x\ — A)5 4 (y) , where the 
vectors x and y span 2D and 4D planes. We have 

A ~ 

1 „ A 1 

O ol „ al = 



JV 2 S xS TRti a ' 

AS « 0.179 xl(T 2 . (B.12) 

A A/ 2 

Intermediate temperatures: 1 <C Ti? <C A -1 / 2 

For temperatures high enough above the deconfinement transition, we have shown 
above that the 9 distribution is gapped. To zeroth order in A it is given by, 

p{9) = 4n(TR) 3 y/(P -P, dl = 2n ^ TK)3 ■ (B.13) 

The intermediate temperature joint eigenvalue densities are obtained from the solu- 
tions listed in the previous subsection by multiplying the eigenvalue densities of the 
previous subsection by 2np(9) and letting the radii A become a function of 9, which 
we will denote r(9), where 

r[e) = 7^kt pW ' (B14) 

where C is now a constant that will depend on the particular solution, as we shall 
see below. 

Combining equations (B.13) and (B.14) we find that the solutions within this 
range of temperature have the ellipsoidal like form: 

^fR r2 + 2 ^ TR ^ 2 = 1 - (B.15) 

The topology of these solutions will no longer be S 6 however, but rather lower di- 
mensional spheres and various singular spaces. 
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The S 1 x S 1 x S 3 solution of the previous section becomes S 1 x S 3 fibred over 
an interval. Topologically, this is a singular space which may be described as 
an S 5 in which a linked S 1 and S 3 have been pinched to a point. The space is 
described by (B.15) together with 

x\ + x\ = \{xl + ... + x l)=r 2 . (B.16) 
cr 

The coefficient in (B.14) is C = 0.557... and as before the ratio a = 1.4344 

The S 1 x S 2 x S 2 solution of the previous section becomes S 2 x S 2 fibred over 
an interval. This is a singular space given by S 5 in which two linked S 2 s are 
pinched to a point. It is described by (B.15) together with 

x\ + x\ + xl = xl + x 2 5 + x 2 6 =r 2 . (B.17) 

The coefficient C = ^(4 - ^2). 

The S 1 x S 1 x S 1 x S 1 solution of the previous section becomes S 1 x S 1 x S 1 
fibred over an interval. It is given by (B.15) together with 

x\ + x\ = x\ + x\ = x\ + x 2 6 = r 2 . (B.18) 

The coefficient C = 0.5628 .... 

The S 4 x S 1 solution becomes a squashed S 5 , given by (B.15) together with 

^ + ... + x2 = r 2 ,x 6 = 0. (B.19) 

The coefficient C = 

The S 1 x S 1 x S 2 solution becomes S 1 x S 2 fibred over an interval. Topologically 
this is an S 4 where a linked S 1 and S 2 have been pinched to a point. It is given 
by (B.15) together with 

x\ + x 2 2 = ^(xl + x\ + xl) = r 2 , x 6 = . (B.20) 

The coefficient C = 0.608 . . . , and as before the ratio a = 1.2404 
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• The S 3 x S 1 solution becomes a squashed S , given by (B.15) together with 

xj + ■ ■■ + x\ = r 2 ,x 5 = x 6 = 0. (B.21) 

The coefficient C=^fi. 

• The S 1 x S 1 x S 1 solution becomes S 1 x S 1 fibred over an interval, which is 
topologically and S 3 with two linked S 1 s pinched to a point. It is given by 
(B.15) together with 

x\ + x\ = x\ + x\ = r 2 , x 5 = x 6 = . (B.22) 
The coefficient C = 0.169 ... . 

• The S 1 x S 2 solution goes over to a squashed S 3 , given by (B.15) together with 

x\ + x\ + x\ = r 2 , X4 = x§ = xq = . (B.23) 
The coefficient C = ^f. 

• The S 1 x S 1 solution goes over to a squashed S 2 , given by (B.15) together with 

x\ + x\ = r 2 , x 3 = ■ ■ ■ = xq = . (B.24) 

The coefficient C = 

High Temperatures: TR ~ A -1 / 2 

In this range of T, recall that the relevant action is 

S = J d9d 6 xp(9,x)(P\x\ 2 + Q9 2 ) 

- \ J dddd'd & xd & x'p(d, x)p(0', x') log (\x - x'\ 2 + (0 - 6') 2 ) . (B.25) 

As the temperature increases the width of the 6 distribution becomes narrower. 
Similarly to the SO (6) symmetric case, we expect a second order phase transition to 
occur at a critical temperature, and for p{9) to collapse to a delta function above that 
temperature. We will now describe R symmetry breaking solutions with 

p{6) = 5(6). (B.26) 
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One finds the solutions and evaluates their actions in the same spirit as in the 
previous subsections. It turns out that the order of the actions evaluated on these 
solutions stay the same and one finds the following ordering: 

(B.27) 

The actions have the following values, putting P = 1 without loss of generality, 
respectively: 0.556775, 0.556853, 0.559185, 0.570093, 0.574554, 0.596574, 0.610025, 
0.653427, 0.846574. We note that the first three solutions in (B.27) are extremely 
close to the value of the maximally symmetric solution S 5 (0.554907). 
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